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ABSTRACT 

It is shown herein that planets with eccentric orbits are more likely to transit than circularly orbiting planets 
with the same semimajor axis by a factor of (1 — e^)~^. If the orbital parameters of discovered transiting planets 
are known, as from follow-up radial velocity observations, then the transit-detected planet population is easily 
debiased of this effect . The duration of a planet 's transit depens upon of its eccentricity and longitude of periastron; 
transits near periastron are shorter, and those near apoastron last longer, for a given impact parameter. If fitting 
for the stellar radius with the other transit parameters, this effect causes a systematic error in the resulting 
measurements. If the stellar radius is instead held fixed at a value measured independently, then it is possible to 
place a lower limit on the planet's eccentricity using photometry alone. Orbital accelerations cause a difference 
in the planet's ingress and egress durations that lead to an asymmetry in the transit lightcurve that could be 
used along with the transit velocity measurement to uniquely measure the planet's eccentricity and longitude of 
periapsis. However, the effect is too small to be measured with current technology. The habitability of transiting 
terrestrial planets found by Kepler depends on those planets' orbital eccentricities. While Kepler will be able to 
place lower limits on those planets' orbital eccentricity, the actual value for any given planet will likely remain 
unknown. 

Subject headings: occultations — planets and satellites: individual (HD80606b, HD147506b) — techniques: photometric 



1. INTRODUCTION 

■ There are presently 21 extrasolar planets known to tran- 
sit their parent stars (http : //exoplanet . eu/). Radial ve- 
locity measurements of all but one of them are consistent 
(with in errors) with circul ar orbits; i.e., zero eccentricity 
{e.g., iLaughlin et al. 2005 ). Presumably, any initial ec- 
fcentricity in th ose orbits has since been damped by tidal 
20001 ). In light of the discovery 



of transiting planets will be found that will not have been 
tidally circularized. B ased on the finding s from radial ve- 
locity planet searches ( Butler et al.l 120061 ). many of these 



circularization (ITrillin 



6f the first transiting extrasolar planet with an eccentric 
orbit, HD147506b, I explore the effect that orbital eccen- 
tricity has on transit lightcurves with an eye toward the 
data to come from CoRo T and Kepler. 

Since tidal circularization is most effective at short 
planet-star distances, as transit search programs extend 
into longer-period regimes the prospects for detecting non- 
circularly orbiting planets grows. Perhaps not coinciden- 
tally, the first known transiting planet on an eccentric or- 
bit, IID147506b (e = 0.507), is t he transiting plane t with 
the longest period (5.63 days) (jBakos et al.l 120071 ). Re- 
cently, Spitzer measurements of the relative timing of the 
secondary eclipse of GJ436b have confirmed that planet's 
nonzero or bital eccentricity an d measured it to be e = 
0.14± .01 (jPemorv et al.ll2007t ). 



newly-discovered transiting planets are likely to follow ec- 
centric orbits. 

Orbital eccentricity has several effects on planetary 
transits. The timing of the transit, relative to that 
of the secondary eclipse, is a strong function of the 
plane t's orbital eccentricity and longitude of periastron 
{e.g., ILaughlin etHI l2005t I Winn et all l2007f ). Planets 
with eccentric orbits, if sufficiently close to their parent 
stars, can have their rotations brought into tidal equilib- 
rium at rotation rates g reater than their mean motions 
("Bar nes fc Fortnev 2003h : this affect s planetar y transit 
lightcurves via the plane t's oblateness (jSeager fc Hui 200^ 
[Barnes fc FortnevI l2003h . These effects and others have 
been explored in the conte xt of eclipsin g binary stars on 
eccentric orbits as well (e.g. iNelson fc D avis 1972). 



With the space-based transi t searches of CoR o T (jBorde et al 



20031 ) and particularly Kepler (jBasri et al.ll2005[ ). hundreds 



In this paper, I explore three additional ways that or- 
bital eccentricity affects the transits of extrasolar planets. 
First, I calculate the increased transit probability for plan- 
ets on eccentric orbits. Next I point out the variability in 
transit duration that results from planets moving faster 
near periastron and slower near apoastron. The third ef- 
fect that I explore is the possibility for asymmetric transit 
lightcurves induced as the planet's trajectory evolves be- 
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tween transit ingress and egress. I then use least-squares 
fits of syntlietic transit liglitcurves to determine whether 
or not these effects can be used to constrain the orbital 
elements of transiting planets from photometry alone. 

Though orbital eccentricity is not inherent to the planet 
itself, I will sometimes refer to planets on eccentric orbits 



as 'ec centric planets' for brevity, following (jLissauer et al 
1993). 



2. EFFECTS 

An extrasolar planet on an eccentric orbit has three pri- 
mary differences relative to that same planet on a circular 
orbit with the same semimajor axis. The eccentric planet 
is more likely to transit, and, if it does, then the transit 
duration depends upon both the impact parameter b and 
the orbital true anomaly /, and the transit lightcurve may 
be asymmetric due to accelerations during the transit. 

2.1. Transit Probability 

Planets on eccentric orbits are more likely to transit 
than equivalent planets with the same semimajor axis but 
circular orbits. Though these planets spend a major- 
ity of their time at greater asterocentric distances than 
their semimajor axes, they spend a majority of their true 
anomalies at smaller asterocentric distances. The proba- 
bili ty for a planet on a c ir cular orbit to transit was derived 
by iBorucki Sz Sumrner i (|1984[) based on the solid angle 
swept out by a planet's shadow: 



P 



(1) 



where p is the transit probability, is the stellar radius, 
and Tp is the distance between the planet and the star. 



Using the method of iBorucki fc Summer i (|l984l ). then, 
the transit probability for an extrasolar planet is equal 
to the solid angle swept out by the planet's shadow, a 
function of both / and the polar angle from the orbit plane 
9, normalized to At: steradians: 



1 
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cos e de df 



(2) 



using 9o and 9i for the angular extent of the shadow below 
and above the orbital plane. Due to the symmetry of the 
problem 9o = — 9i. Geometry sets 9i = sin~^ (see 
Figure [1]). Hence, integrating over 9, 
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and plugging in and 9i 



sin 6* I ' df 
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(3) 



(4) 



For planets on eccentric orbits, rp varies with time. 
However, the variation of function of / is all that 



matters for determining the solid angle over which the 
planet will transit. 



ap(l-e2) 



(5) 



1 + e cos / 

( Murrav fc Dermott where e is the planet's orbital 

eccentricity and Op its semimajor axis. 

Plugging in rp from Equation [5] and integrating over / 
2i?* (l + e cos/) 
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d/. 



(6) 



leads to 



2i?* 



P ^ 1- 



Att ap(l - e2) 



2tt 



l + ecos(/) df. 



(7) 



which leads to the result that for planets on eccentric orbits 



ap(l-e2) 



(8) 



since f^^ df = 2tt and J^^ cos f df = 0. 



The above derivation is valid for transit impact param- 
eters 6 <= 1. To exclude all grazing transits, replace the 
i?* in the numerator of Equation [5] with i?* — Rp, where 
Rp is the radius of the transiting planet. Similarly, to in- 
clude all transits, no matter how grazing, the numerator 
of Equation [8] would be i?, + Rp. 

The increased transit probability for even a planet with 
a significantly e ccentric orbit with e = 0.5, similar to that 
for HD147506b (jBakos et al.ll2007l ). is modest: 33%. How- 
ever, the increased probability fo r planets on extr emely 



eccentric orbits hke HD80606b (|Naef et al.l l200lh with 



e = 0.93 is 640%! Given that 28 of the 224 planets with ra- 
dial velocity orbits have e > 0.5 (http : / / exoplanet . eu/[ ), 
transit surveys should detect a significant number of plan- 
ets on eccentric orbits. Nearly half, 110 out of 224, of radial 
velocity planets are more eccentric than our solar system's 
most eccentric planet. Mercury (e = 0.2056). Planets with 
extreme orbital eccentricities will be detected at a rate de- 
cidedly higher than their occurrance would predict given 
Equation [TJ 

This excess will lead to a bias in the raw planet inci- 
dence as a function of semimajor axis derived from CoRoT 
and Kepler discoveries. The bias can be easily corrected 
by accounting for the (1 — e^)~^ detection increase factor 
from Equation [51 or , more formally, by using the original 
IBorucki fc Summer^ ( 19841 ) probability (Equation[T]) while 
substituting the instantaneous planet-star distance at mid- 
transit for the orbital radius for each detected planet. 

Debiasing requires knowledge of the planet's orbital ec- 
centricity and longitude of periapsis, which would be diffi- 
cult to ascertain for planets too small to induce detectable 
radial velocity variations (Section [3]). Although this ex- 
tra step adds complexity, it is heartening to note that the 
orbital eccentricity bias induced in planet distrubutions as 
derived using the transit method is precisely calculable and 
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removable. The orbital eccentricity bias in radial velocity 
planet surveys (as results from data gaps near periastron 
passage, for instance) is known less precisely and is there- 
fore more challenging to remove. 

The differing transit probability at periastron and 
apoastron can lead to planets that transit, but have no 
secondary eclipse if the planet is sufficiently eccentric, in- 
clined, and transits near periastron. Alternately there can 
exist planets that show secondary eclipses but no primary 
transit. In the case where giant planets have just sec- 
ondary eclipses, the small secondary eclipse depth could 
be mistaken for the primary transit of a terrestrial-sized 
planet. Careful monit oring of reflected- light phase vari- 
ability (as described in Jenkins! ( 20021 )) may eliminate this 
source of systematic error. Secondary-eclipse-only planets 
would need to be very near their parent stars during sec- 
ondary eclipse in order to have sufficient reflected light so 
as to mimic the transit of a terrestrial-sized planet. Faint 
secondary stars in similar orbits could mimic the terres- 
trial planet transit regardless of their distance from the 
primary star at secondary eclipse. 

2.2. Transit Duration 

For a given star mass, planets with the same orbital 
semimajor axes have the same energy per unit mass 
{—GMtf/2ap)] those with higher eccentricities have lower 
specific angular momenta. As such, those more highly 
eccentric objects move faster near periastron, and slower 
near apoastron. If such an eccentric planet were to transit, 
it would then have transits of shorter or longer duration 
than the equivalent planet on a circular orbit transiting 
with the same impact parameter (6). 

The velocity of a circularly orbiting planet Vcirc is 



(9) 



Using conservation of energy, the periapsis velocity Vpori 
of a planet with orbital eccentricity e can be shown to be 
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cri 



pcricirc 



(10) 

In the case of planets discovered by their transits, Op is 
set by the planet's known orbital period, and hence the 
comparison to Vdrc is most relevant. For convenience I also 



compare Fperi to Vp 



the velocity of a planet orbiting 



circularly at the eccentric planet's periastron. 
at apoapsis 



Similarly, 



1 



1 



for similarly named variables (see Figure [5)) ■ 



yl ^ ^^pocirc 

(11) 



The pericenter and apocenter velocities behave as ex- 
pected in their extremes. Vpeii approaches escape velocity 




Fig. 1. — Geometry of a planetary orbit and the 
solid angle swept o ut by its shadow, modified from 
Borucki fc Summer^ ( 1984 ) to account for orbital eccen- 
tricity. 




Fig. 2. — This figure illustrates the variables described in 
Section [2T2I of the text, as applied to a hypothetical planet 
with e = 0.5. 
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(V2 Fpe 



as e — !■ 1 but is undefined for non-closed or- 



bits where e > 1. Va 



apo 



as e — > 1, in absolute terms, 



as a function of Vcirc, and as a function of Vapocirc- 

A significantly eccentric planet with e = 0.5 (similar 
to HD147506b for which e = 0.507) travels (- 1.73) 
times faster at periapsis than you would expect given its 
period and assuming a circular orbit, and times more 
slowly at apoapsis. This planet's transit duration, if it 
were to transit at periastron, would be 58% as long for a 
particular impact parameter as the transit of its circularly 
orbiting equivalent. If transiting at apoastron (three times 
less probable than a periastron transit; see Section [27T|) . the 
transit would last 73% longer than the equivalent circularly 
orbiting planet. 

A plane t on an extremely eccentric orbit like HD80606b 
(e = 0.93) (jNaef et al.ll200ir ) would have a periastron tran- 
sit duration only 19% that of its same-semimajor-axis cir- 
cular equivalent, and 72% that of a planet orbiting circu- 
larly at HD80606b's periastron. Conversely if HD80606b 
were to transit near apoastron, such a transit would last 
5.25 times longer than if HD80606b were to transit in a cir- 
cular orbit with the same semimajor axis, and 3.78 times 
longer than if HD80606b were to transit in a circular orbit 
at its true apoastron. 

The duration of a cxoplanetary transit depends on the 
chord length of the planet's apparent passage in front of 
the star (2i?«cos6) and the planet's azimuthal velocity 
Vf = rf. In the cases discussed above, when a planet is 
at periastron and apoastron, Vpci-i and Vapo are equal to 
Vf. The rest of the time, Vf if not equal to the planet's full 
velocity as the planet will also have a velocity component 
radial to the star (K)- According to iMurrav fc Dermott 
(|2000l) . Vf varies sinusoidally with /: 



Vf 



Kir 



1 + e cos / 



(12) 



Hence the planet's velocity is greater than the equivalent 
circular orbit velocity for more than half of the orbit as a 
function of the true anomaly. 

2.3. Transit Symmetry 

The dependence of the planet's azimuthal velocity on / 
as shown in Equation[T2]belies to another effect that eccen- 
tric orbits have on transits. Because Vf changes slightly be- 
tween ingress and egress (unless the planet is at periastron 
or apoastron mid-transit), eccentric planets can produce 
asymmetric transit lightcurves. If a planet transits after 
periastron and before apoastron, the time that it takes for 
the planet to ingress across the star's limb will be shorter 
than the time that it takes to egress. Similarly if the planet 
is between apoastron and periastron, then the ingress will 
be longer than the egress. 

To calculate the velocity difference between ingress and 
egress, I first take Vf from Equation [12] and differentiate 



it with respect to /: 

dVf_ 
df 



e Vc w 



sin/. 



(13) 



The total ingress-egress velocity difference, AV, is equal to 
times A/, the difference of the true anomaly between 

ingress and egress, under the assumption that varies 
negligibly across the transit. Taking /q to be the true 
anomaly of the planet at mid-transit, 



A/ = 



2i?, 



(14) 



assuming that Tp varies only slowly during the transit. 
Now, plugging in Tp from Equation O 



AV^ = ^ A/ 



2R^eVc 



-sm/o (1 + ecos/o) 



a(l - e2)3/2 



(15) 

I show a plot of the varying component of AV for various 
values of e in Figure[31 The fractional difference in velocity, 
AV/V, is a bit simpler: 



AV 
V 



2R^e 
a(l-e2) 



sin/o , 



(16) 



and has an evident maximum where /q = ±Tr/2. Hence, 
the greatest fractional variation in ingress and egress du- 
ration will occur ninety degrees away from periastron and 
apoastron, as measured in the planet's true anomaly. 

To determine where the absolute velocity difference AV 
is maximized, I differentiate with respect to / and set 
the result equal to zero: 







dAV 



2i?,eVcirc d [ sin /o -he sin/o cos /o] 



a(l — e^) 2 



(17) 

Differentiating and using the double-angle formula sin(2x) = 
2 sin X cos a;, I determine that 



cos/o + ecos(2/o) 







(18) 



I show the solution to this equation in Figure 3) The max- 
imum velocity difference occurs at /o — ±7r/2 for infinites- 
simal eccentricities, and approaches /o = ±7r/3 as e ^ 1. 

To illustrate the consequences of this effect, I apply the 
results derived above to newly-discovered e ccentric tran- 
siting planet HD147506b (iBakos et al.ll2007l). HD1475 06b 



was determined by its discoverers. iBakos et al.l (l2007|), to 



have a radius of 1.18 Rjup, an orbital semimajor axis of 
0.0685 AU, an orbital eccentricity of 0.507, and a longi- 
tude of periastron of 184.6°. The planet's parent star, 
HD147506, was determined to have a radius of 1.8 Rq and 
a mass of 1.35 M© (note that refined system parameters 
were published by ? while this paper was in review - the 
new values do not change the qualitative results that I 
describe here, but future work should employ these newer 
values) . I show a to-scale schematic of the system in Figure 
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Fig. 3. — Velocity difference between egress and ingress 
as a function of the planet's orbital true anomaly at mid- 
transit, fo- I have only plotted that portion of Equation 
[15] that is a function of / and e; to convert to true AV, 
multiply by 2E*ysin^_ As an example, the multiplier for a 
planet in an Cp = O.lAU orbit around a IMq star is 8.76 
km/s. 
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Fig. 4. — The true anomaly at mid-transit (/o, here in 
units of radians divided by tt) for which the difference be- 
tween a planet's transit ingress and egress velocity is max- 
imized as a function of the planet's orbital eccentricity, e. 



From Equation [12] HD147506b's Vcirc is 132.2 km s^^. 
This planet's orbital parameters are particularly favor- 
able with respect to maximizing the magnitude of A Vf. 
HD147506b's large eccentricity, fast Kirc, and nearly ideal 
true anomaly at midtransit (/o — 4.84, very close to 
— 7r/2), combined with HD147506's large stellar radius 
yield A Vf = — 24.6 km s^^. Calculating the ingress 
and egress times r using 

r = (19) 

Vf cos(sm b) 

with the appropriate Vf values for the transit ingress and 



egress of HD147506b yields Ti, 
T„„,.„=, = 8 min 37 sec. 



DETECTABILITY 



10 min 36 sec and 



The CoRoT and Kepler missions will d iscover hun- 



dreds of new transitin g extrasolar planets (jBorde et al 



2003 : Basri et al. 2005f l. The most massive of these will 



be amenable to radial velocity follow-up observations to 
measure their masses; a more complete set of radial veloc- 
ity measurements (i.e., covering the full orbit and not just 
the radial velocity maxima and minima, which are all that 
would be required to determine mass assuming a circular 
orbit and the epoch and period as established by the tran- 
sit) can determine the planets' orbital eccentricities and 
longitudes of periastron. However for planets of Neptune- 
mass and smaller, in larger orbits, or orbiting fainter stars, 
radial velocity follow-up will not be practical or in some 
cases possible. In those cases, it would be useful to at- 
tempt to constrain the orbital eccentricity of those plan- 
ets using transit photometry alone. Such a determination 
would bear critically on the climatic variability of Earth- 
like worlds; highly eccentric planets may not be habitable 
even if their orbital semimajor axes place them within their 
stars' habitable zone. 

To determine the photometric detectability of the tran- 
sit duration and asymmetry effects of eccentric planet or- 
bits, I create synthetic transit lightcurves that I then fit 
as if the orbit were circular. I assume knowledge of the 
planet's period, as that value will be measured by the 
time between transits in Kepler and CoRoT data. I cal- 
culate both the synthetic lightcurves and the best-fit solu- 
tions using th e analytical approxima tion of lMandel &: Ago! 
([2002). The iMandel fc Agoll (|2002l ) formulation assumes 
that the portion of the star covered by the planet has uni- 
form surface brightness, but that brightness accounts for 
limb-darkening; hence it is least accurate for ingress and 
egress. However, the detrimental effects on this particu- 
lar calculation should be minimal since both the synthetic 
and best-fit lightcurves should show the same systematic 
errors, which, when subtracted, should leave a good esti- 
mate for the proper fit res i duals . The effects of light-travel 
delay illustrated bv iLoebl (|2005l ) are not included. 



For larger planets, those where the ingress and egress 
are temporally resolved, that have high signal-to-noise 
lightcurves, I first fit for R^^ Rp, the transit impact pa- 
rameter b, and a stellar limb darkening coefficient ci, the 
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treatment applied to HD209458b by iBrown et all (|200ll ). 
I assumed a 1 i?jup planet orbiting a 1 Rq, 1 Mq star 
with Op = 0.1 AU and e = 0.5. Since Rp/R* is set by the 
transit depth, and b by the ingress/egress time relative to 
the total transit duration, orbital eccentricity in this type 
of fit systematically affects i?, , which governs the total 
transit timescale. In my test runs, the best-fit i?*moasurod 



varies as ii* 



^measured = i?*VcircV^(, • Hence assuming a cir- 
cular orbit when the orbit is actually eccentric leads to a 
systematic error in the measurement of i?* and, by exten- 
sion, Rp [Rp/Ri, is unaffected). The measured values are 
smaller than the actual values if the planet is near perias- 
tron, and are larger if the planet is near apoastron. 

This systematic error can be addressed by assuming a 
stellar radius, as could be estimated by other observations 
such as parallax, spectral type, and stellar absolute mag- 
nitude. In this case, the transit timescale can be set by 
fitting explicitly for Vf^ in addition to rp, 5, and ci. How- 
ever, without knowledge of the mid-transit true anomaly 
/o, Vfo cannot uniquely determine the orbital eccentricity. 
Instead, Vfg can constrain e if we allow that the planet 
must have a minimum eccentricity in order that Vf reach 
IfV/o > K;irc, then from Equation [To] 



e > 



^ Vcirc 



1 



(20) 



and if Vjg < Vciic , then similarly from Equation [TT] 

1 



e > 



1 



- Voire 



(21) 
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Fig. 5. — To-scale diagram of the HD147506 system. 
Ingress and egress times listed are in minutes : seconds 
format. Figure concept inspired by Gregory Laughlin's 
[http : //oklo ■ org/, website. 



The above lower limits can only be placed when the 
planet's ingress and egress are resolved. In order to con- 
strain the eccentricity of terrestrial-sized planets, both 
high temporal cadence and high photometric precision 
would be necessary. Resolving the ingress of the Earth, 
with Tingicss = 7.04 minutes at 6 = 0, will not be possi- 
ble given Kepler data alone. Coaddition of multiple tran- 
sits from multiple telescopes might help to constrain the 
ingress times for detected transiting terrestrial planets. 

The degeneracy between e and /o, multiple comina- 
tions of which can produce the same Vjg , can be broken 
by measuring the transit asymmetry outlined in Section 
12.31 To measure the detectability of the asymmetry, I pro- 
duce a synthetic lightcurve for HD147506b (rp — l.lSRjup, 
i?, = 1.8Rq, an = 0.0685 AU. e = 0.507, /o = 4.84, 
and assuming iBrown et al. ( 2001 ) limb darkening coeffi- 
cient ci = 0.64) that I fit using a model planet with a 
circular orbit. The residuals, which I refer to as the de- 
tectability, are shown in Figure [H] for transits at several 
impact parameters. 

When faced with an ingress and egress of differing dura- 
tion, the best-fit circular-orbit planet model splits the dif- 
ference. On ingress for HD147506b, which transits before 
periastron (/q = —1.44), the real planet is moving more 
slowly than the model planet. In the difference lightcurve, 
the real planet hits the stellar limb before the model. 
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Fig. 6. — Detectability of the transit lightcurve asymmetry 
induced by the orbital eccentricity of planet HD147506b, 
for impact parameters b = 0.0, b = 0.3, b = 0.7, and 
b = 0.85. 
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such that at mid-ingress the real and model planets are 
in the same place. So the real-minus-model detectability 
initially trends negative, reaches zero at mid-ingress, and 
then trends positive until the end of the planet's ingress. 
The detectability varies slowly between the end of ingress 
and the beginning of egress, resulting from the planet cov- 
ering stellar areas that have been slightly differently limb- 
darkened. 

The process is reversed on egress. With the planet 
now moving more quickly than the circular-orbit model, 
the model begins its egress early so that the real planet 
will have caught up with it by mid-egress. Hence the 
real-minus-model detectability is again initially negative, 
reaches zero by mid-egress, and then trends positive until 
the end of egress for both the real and model planets. 

Despite the large AV for HD147506b, the magnitude of 
the transit lightcurve asymmetry induced by orbital eccen- 
tricity is rather low (Figure [6]), peaking at only 3 x 10~^ 
of the stellar flux. This low detectability is unmeasurable 
given current capabilities, and will probably prove chal- 
lenging in the future as well given inherent stellar vari- 
ability. The large stellar radius of HD147506 increases 
AV, but also decreases i?p/i?*, diminishing the asymme- 
try effect. Hence though measurement of the asymmetry 
induced by an eccentric planet orbit can, along with Vf, 
uniquely determine both e and the longitude of periastron, 
the small magnitude and duration of the effect are such 
that a measurement is unlikely to be practical. 

For transiting terrestrial-sized planets, determination of 
orbital eccentricity is even more difficult. If the ingress and 
egress are temporally resolved, that duration along with 
the total transit duration together set the transit impact 
parameter and Vf, providing a minimum e constraint as 
per Equations [20l and [2T] For objects with such small tran- 
sit depths as terrestrial planets, though, it will be difficult 
to temporally resolve the planets' ingress and egress. 

Another way to constrain b for terrestrial planets 
would be to use the ef f ect of stellar limb darkening. 
Seager fc Mallen-OrnelasI ( 20031 ) showed that minimiz- 
ing limb darkening more precisely delineates the end of 
a planet's ingress and the beginning of its egress. However 
strong limb darkening, if well-understood, can provide a 
mechanism to ascertain a tr ansit's impact parame ter. Sim- 



ilar to the analysis used by iKnutson et al.l (|2007r ) to study 
HD209458b, fixing limb darkening coefficients based on 
theoretical or previousely determined values could suffice 
to constrain b, which, with a previously measured i?*, 
would then determine Vf and allow for constraints on e. 

4. CONCLUSIONS 

For a given orbital semimajor axis, extrasolar planets 
on eccentric orbits are more likely to transit than plan- 
ets on circular orbits by a factor of (1 — e^)^^. As the 
space-based transit surveys CoRoT and Kepler discover 
transiting planets that are far enough from their parent 
stars to have avoided tidal circularization, more highly ec- 



centric planets will be found preferentially. This bias is 
straightforward to remove if the eccentricity and longitude 
of periastron are known, as they could be given follow-up 
radial velocity observations. 

The duration of a transit is a function of the planet's 
tangential velocity at mid-transit, Vf^. For eccentric plan- 
ets Vfg is greatest at periastron and smallest at apoastron. 
Hence if a planet transits near periastron the duration is 
shorter than that of an equivalent planet in a circular or- 
bit, and similarly transits that occur when a planet is near 
apoastron is longer than those of the equivalent circularly 
orbiting planet. It would be useful to be able to use this 
effect to determine the orbital eccentricity of discovered 
transiting extrasolar planets, either before or without ra- 
dial velocity follow-up. 

If fitting the resulting lightcurve with a model planet 
on a circular orbit with the known planetary period, a 
systematic error results in the determination of the tran- 
sit paramete r s if fi tting for i?,, i?p, b, and ci as done by 



Brown et al.l (|200ll ) for HD209458b. If instead the model 
system uses an assumed stellar radius measured by differ- 
ent means, then the transit velocity Vf^ can be measured. 
However without another, independent measurement of ei- 
ther e or the planet's longitude of periastron, knowledge of 
Vfg cannot alone determine those parameters. It can set a 
lower limit on a planet's orbital eccentricity. 

The difference in Vf between a planet's ingress and 
egress that results from the planet's orbital accelerations 
can resolve the e / longitude of periastron degeneracy. 
This velocity differential AV introduces an asymmetry into 
the transit lightcurve: ingress is longer than egress be- 
fore periastron, and shorter after periastron. My model 
fits to synthetic eccentric planet transit lightcurves show 
that the detectability of this asymmetry is small, of order 
3 X 10~^ for recently discovered eccentric transiting planet 
HD147506b. An effect that small is undetectable using 
present techniques. As HD147506b is nearly a model can- 
didate for which to observe this effect, it is unlikely that 
transit lightcurve asymmetry will prove useful for deter- 
mining the orbital parameters of transiting planets using 
photometry alone. 

Determination of orbital eccentricity is of critical im- 
portance for evaluating the habitability of terrestrial-sized 
transiting planets discovered by Kepler. As these planets 
have masses too low for radial velocity measurements to 
detect, our only constraints on e for these planets will be 
photometric. If the stellar radius can be assumed from 
other, prior measurements, then it is possible to use theo- 
retical stellar limb-darkening coefficients within the Kepler 
bandpass to measure the transit impact parameter. This 
measurement would then constrain V^,, and allow a lower 
limit to be placed on e. 

No techniques currently available will be able to 
uniquely measure the orbital eccentricity of the terres- 
trial extrasolar planets that Kepler will discover. The 
lower limits on eccentricity described above will allow for 
a statistical exploration of the eccentricity distibution of 
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terrestrial planets. That distribution will serve to con- 
strain the formation and evolution of such planets, as it 
has done for giant planets. However, whether or not any 
particular Kepler planet is truly habitable will remain 
unknown until its orbital eccentricity can be measured. 

JWB acknowledges the support of the NASA Postdoc- 
toral Program, administered for NASA by Oak Ridge As- 
sociated Universities, and the support of NASA's Kepler 
mission for publication costs. 
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